Glassy dynamics of partially pinned fluids: an alternative mode-coupling approach 
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We use a simple mode-coupling approach to investigate glassy dynamics of partially pinned fluid 
systems. Our approach is different from the mode-coupling theory developed by Krakoviack [Phys. 
Rev. Lett. 94, 065703 (2005), Phys. Rev. E 84, 050501(R) (2011)]. In contrast to Krakoviack's 
theory, our approach predicts a random pinning glass transition scenario that is qualitatively the 
same as the scenario obtained using a mean-field analysis of the spherical p-spin model and a mean- 
field version of the random first-order transition theory. We use our approach to calculate quantities 
which are often considered to be indicators of growing dynamic correlations and static point-to-set 
correlations. We find that the so-called static overlap is dominated by the simple, low pinning 
fraction contribution. Thus, at least for randomly pinned fluid systems, only a careful quantitative 
analysis of simulation results can reveal genuine, many-body point-to-set correlations. 



PACS numbers: 64.70.P, 64.70.Q, 61.20.Lc, 05.20Jj 

Recently there have been several theoretical and sim- 
ulational studies of glassy dynamics of fluid systems in 
which some particles, randomly selected out of an equilib- 
rium configuration, have been frozen or pinned [TH51 [5] . 
Originally these so-called partially pinned systems were 
considered to be just one special example of a broad class 
of model porous systems known as quenched-annealed bi- 
nary mixtures [5141 lj. However, it has now been realized 
that glassy partially pinned systems can be used to reveal 
still unresolved aspects of the glass transition. 

First, it was proposed that by analyzing systems in 
which some particles, taken out of an equilibrium con- 
figuration, have been frozen, one can study a growing 
"amorphous order" that is supposed to develop in glassy 
fluids |12j . In early studies this idea was implemented 
using the so-called cavity geometry: all particles except 
those within a spherical cavity were frozen and the lo- 
cal overlap of the original equilibrium configuration with 
configurations equilibrated in the presence of pinned par- 
ticles was monitored [T3]. It was argued that the de- 
pendence of this overlap on the cavity diameter reveals 
a length characterizing the so-called static point-to-set 
correlations, i.e. correlations between the density at the 
center of the cavity (point) and the positions of the frozen 
particles (set). It was shown [T3] that, at least in simple 
models, these point-to-set correlations grow with increas- 
ing relaxation time. Subsequently, other geometries were 
introduced: one in which all particles except ones in a 
layer are frozen (the sandwich geometry), one in which 
all particles in a semi-infinite space are frozen (the wall 
geometry) |15j or one in which a randomly chosen subset 
of particles, distributed uniformly throughout the system 
is frozen. It has been argued [5] that the last geometry, 
i. e. the partially pinned system, is the best candidate to 
study growing static correlations. 

The second motivation for the recent interest in par- 
tially pinned systems comes from the realization that 
pinned particles, while maintaining the equilibrium struc- 



ture of the fluid [10l COS] , may induce an ideal glass transi- 
tion at temperatures or densities that are more accessible 
to computer simulation studies [3] . Thus, the analysis of 
this so-called random pinning glass transition 8\ could 
both shed light on the glass transition itself and provide 
a new way to test diverse theoretical descriptions used to 
describe it. In particular, Krakoviack j4] recently argued 
that two different approaches, the mode-coupling theory 
[T7] and the random first order theory (RFOT) 18! make 
strikingly different predictions for the glassy behavior of 
partially pinned systems. Thus, a simulational study of 
glassy behavior of a partially pinned system could easily 
disprove one of these approaches. 

The conclusion reached in Ref. [3] was surprising for 
two reasons. First, it is usually assumed that mode- 
coupling and RFOT theories are related. In particular, 
the former theory is supposed to describe rather well the 
onset of glassy behavior while the RFOT approach is 
supposed to replace the mode-coupling transition with 
a crossover and to describe properties of deeply super- 
cooled fluids. It has to be admitted that recent theo- 
retical analyzes of these theories in higher dimensions 
revealed some rather disturbing discrepancies [19H21j . In 
spite of this fact, qualitative disagreement between them 
in three dimensions was not expected. Second, Cam- 
marota and Biroli [5] showed that a mean-field analysis 
of the p-spin model with partially pinned spins predicts 
qualitatively the same random pinning glass transition 
scenario as the mean-field version of the RFOT theory. 
Thus, we are now faced with a rather unpleasant qual- 
itative disagreement between two mean-field-like calcu- 
lations, mode-coupling theory and the mean-field p-spin 
model. This disagreement is even more striking if we 
recall that (in the case of an un-pinned system) the so- 
called schematic model of mode-coupling theory is iden- 
tical to the p-spin model (for p = 3 used in Ref. [5]). 

Faced with the above described conundrum, we shall 
recall that there is some freedom in the formulation of 
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the mode-coupling approach, especially for more com- 
plex systems like mixtures and partially pinned systems. 
Our goal in this Letter is to propose a simple, alternative 
mode-coupling approach that results in the random pin- 
ning glass transition scenario which is qualitatively con- 
sistent with the mean-field analysis of both the p-spin 
model [5] and the RFOT theory [3]. In addition, we will 
use our approach to calculate quantities that have been 
used in earlier simulational studies to monitor the growth 
of dynamic correlations and static point-to-set correla- 
tions. We will show that a careful analysis of these quan- 
tities is required to reveal genuine many-body effects. 

A mode-coupling theory for a partially pinned system 
can be derived in different ways. Here, we will outline a 
projection operator derivation which is easily compared 
with Krakoviack's theory [31 [52]. Also, following Refs. 
[31 [22] we will refer to the mobile (un-pinned) particles 
as the fluid particles and to the pinned ones as the ma- 
trix particles. The fundamental dynamical variable used 
in our theory is the Fourier transform of the microscopic 
fluid density, n-^(q) = Y^jli e~^' rj , where Tj denotes the 
position of the jth fluid particle and Nf is the number 
of fluid particles. To describe the dynamics of the sys- 
tem, we use the fluid intermediate scattering function, 
F(q;t) — j^j (n-^(q) exp(Qt)n f (— q)), with being the 
system's evolution operator and (...) denoting the average 
over all (fluid and matrix) particles of the system. For 
simplicity, we assume here that the microscopic dynam- 
ics is Brownian, thus f2 is the many-body Smoluchowski 
operator. To derive an equation of motion for F(q; t) 
we follow the standard procedure (T71 [23] and arrive at 
an exact but formal memory function representation for 
F(q; t) [23]. We project the so-called fluctuating force on 
the space spanned by the products of the fluid densities, 

(k)nt (q — k), and the products of the fluid and matrix 
densities, r?/(k)n m (q— k), where n m (q) = Ylj=i e _lq Sj 
is the density of the matrix (pinned) particles, with s-,- be- 
ing the position of the jth matrix particle and N rn being 
the number of the matrix particles. Finally, we factorize 
the resulting four-point correlation functions [25 . 

We obtain the following equation of motion for F(q; t), 

du (6{t -u) + M ilT (q; t - u)) d u F(q; u) = 

Dq 2 n f n m c(q)h{q) - Dq 2 {l - n f c(q))F(q;t). (I) 

Here D is the diffusion coefficient of an isolated fluid 
particle, n f and n m are the densities of the fluid and ma- 
trix particles, respectively. h(q) and c(q) are the Fourier 
transforms of the complete system's correlation function 
and direct correlation function, respectively |26j . Finally, 
M lrr (g;t) is the irreducible memory function, 



where 



2 >(q, k) = ^ [q • kc(fc) + q • (q - k)c(|q - k|)] 2 (3) 

V«(q, k) = n m D [q • (q - k)c(|q - k|)] 2 

x [I + n m h(\q - k|)] + 2n f n m Dq ■ (q - k)c(|q - k|) 

x [q • kc(fc) + q • (q - k)c(|q - k|)] fc(|q - k|) (4) 

V (0) (q,k) = Z>n^q-kc(fc)q-(q-k)c(|q-k|) 

x(l-a;)/i(fc)ft(|q-k|) (5) 

In Eq. @, x is the fraction of particles that are pinned, 
x = n m /n, with n = nj + n m being the total density 
of the system (in Refs. [SI El IB] the pinning fraction is 
denoted by c; we changed the notation to avoid confusing 
the pinning fraction and the direct correlation function) . 

It should be noted that there is a time-independent 
term at the right-hand-side of Eq. ([!]) and there is also a 
time-independent contribution to M lrr originating from 
V~(°\ Eq. |5]). These two terms make the long-time limit 
of F(q; t) non-zero even below the glass transition. 

Before turning to the results obtained from Eqs. ^ 
[5) we compare our present approach to Krakoviack's 
[HI22]- Krakoviack also follows the projection operator 
procedure detailed in Ref. [17] . However, he takes as 
the fundamental dynamical variable the so-called relax- 
ing part of the fluid density, <5n/(q) = n-^(q) — (n(q))^ 
where (...) denotes the equilibrium average over the 
positions of the fluid particles only, with the matrix 
(pinned) particles treated as a set of fixed obstacles. Fur- 
thermore, he projects the fluctuating force on the space 
spanned by the products of the relaxing fluid densities, 
Sri? (k)5nf (q— k), the products of the relaxing fluid den- 
sity and the matrix density, <5n^(k)n m (q — k), and the 
products of the relaxing fluid density and the average 
fluid density for a given set of positions of the matrix par- 
ticles, Jn-^k) (n(q — k)) . It should be noted that, while 
the all the variables that we use are either one-particle 
additive or pairwise-additive, the average fluid density 
(n(q))^ used in Ref. [52] is not one-particle or pairwise 
additive in terms of the matrix particles; instead, it in- 
cludes terms involving arbitrarily many matrix particles. 

The main difference between our approach and that 
of Refs. [3] H2] is that we use the mode-coupling ap- 
proach to find both the time evolution of the relaxing 
density fluctuations and the average non-relaxing den- 
sity fluctuations. In our language the latter quantity is 
the long-time limit of the fluid intermediate scattering 
function, F(q) = lim^oo F(q; t). In contrast, Krakovi- 
ack uses the mode-coupling approach only to find the 
time evolution of the relaxing density fluctuations and 
resorts to a separate static approach, replica OZ inte- 
gral equation theory [SI [TU], to calculate the proper- 
ties of the non-relaxing density fluctuations. In his lan- 
guage the latter fluctuations are characterized by the 
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+W 1 )(q,k)F(fc;t) + y(°)(q,k)l (2) 
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so-called disconnected or blocked fluid structure factor, 
S b (q) — jjj(n(q))^ (n{— q))^ where 7TT denotes the aver- 
age over the positions of the pinned particles. 

Two arguments were put forward for the approach used 
in Ref. [52]. First, it was argued that since the calcula- 
tion of the properties of non-relaxing fluid density fluc- 
tuations is a static problem, it is appropriate to use the 
mode-coupling approach for the dynamics of the relax- 
ing part only. Second, in the formulation of Ref. [22] 
the so-called connected fluid structure factor, S c (q) = 
j^j (<5n^(q)5n^(— q)), appears naturally in the expres- 
sion for the so-called characteristic frequency [22] which, 
for a Brownian system, is given by 0J 2 (q) = Dq 2 /S c (q). 

To answer these arguments we note that the mode- 
coupling calculation of the non-decaying part of the fluid 
intermediate scattering function F(q;t) can be given a 
static interpretation [27] . Thus, our approach can be 
considered to be a combination of a specific static cal- 
culation that is consistent with the mode-coupling pre- 
diction for the non-decaying part of F (g; t) and the 
mode-coupling calculation of the decaying part of F(q; t). 
Moreover, if the equation of motion is re- written in terms 
of the decaying part only, the characteristic frequency 
acquires the form w 2 (g) = Dq 2 / '(S-' (g) — F(q)), where 
S (q) = F(q; t = 0) = l + rifh(q) is the fluid static struc- 
ture factor. This form, upon identification of F(q; oo) 
and S b (q), is identical to that obtained in Ref. [22]- 

Eqs. (I][5) allow us to calculate both F(q), i.e. the 
non-decaying long time limit of F(q; t), and the time de- 
pendence of F(q; t). In fact, for the fluid states the small 
x limit of F(q) can be easily obtained from Eq. 



F(q) = xn 2 h 2 (q) + 0(x 2 ). 



It can be showed that this result coincides with the exact 
small x limit of the blocking structure factor, S h (q). 

To solve Eqs. (I~f5l numerically we disctretized the 
wave- vector space using 500 wave- vectors with the small- 
est one equal to 0.05 and a spacing of 0.1. The resulting 
set of coupled integro-differential equations was solved 
using the procedure outlined in Ref. 28 . We used the 
hard sphere interaction potential and the Percus-Yevick 
approximation for the static correlations, h{q) and c(g). 
The results below are presented using reduced units: dis- 
tance is measured in terms of the hard sphere diameter 
a and time in terms of a 2 /D. 

The random pinning glass transition scenario is pre- 
sented in Figs, [IJ and [2] The standard mode-coupling 
transition present at x — extends into the x > plane. 
However, the discontinuity of F(q) at this transition de- 
creases with increasing x and the transition disappears 
at x c rj 0.1395. Beyond x c , the intermediate scattering 
function changes continuously with the volume fraction 
ip = n7r/6. This scenario is qualitatively consistent with 
the p-spin model results presented in Ref. [8] and with 
the mean-filed RFOT analysis presented in ;3L Thus, 
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FIG. 1: (a) dynamic phase diagram of a partially pinned 
hard sphere system, x and <p denote, respectively, the frac- 
tion of particles that are pinned and the volume fraction of 
all (mobile and pinned) particles. Solid line: discontinuous 
mode-coupling transition. Dashed line: localization transi- 
tion (i.e. vanishing self-diffusion coefficient) 29 . The filled 
circle and the thin vertical line indicate states for which the 
time dependence of the intermediate scattering function is 
shown in Fig. [5] (b) Circles: x dependence of the discon- 
tinuity of the long-time limit of the normalized intermediate 
scattering function, f(q) = F(q)/F(q;t = 0), along the mode- 
coupling transition line shown in (a); thin solid line shows the 
square root singularity, Af oc (x — x c ) 1//2 . 
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FIG. 2: Time dependence of the normalized intermediate 
scattering function, f(q;t) = F(q;t)/F(q;t = 0). Thick solid 
lines: x — 0.05 and <p near the mode-coupling transition, solid 
circle in the left panel of Fig. [l] Thin solid lines: x — 0.15 
and ip along the thin vertical line in Fig. [l^,. 



it agrees with the striking prediction of Cammarota and 
Biroli that it is possible to enter the glass phase without 
ever encountering a divergence of the relaxation time. 

We used the numerical solution of Eqs. (l|5l to cal- 
culate the so-called collective and single-particle overlap 
functions. To define these functions we followed Berthier 
and Kob [5]. Specifically, Q(t) is proportional to the the 
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probability that a spherical cell of radius a is occupied at 
both t = and t —> oo, 

Q(*) = 4^3 / dr i / ^(Irx-raht^o-riJ^o-ra) 

+Qrand (7) 

where F(r;t) is the inverse Fourier transform of F(q;t) 
and Q r and = 47ma 3 /3. In the definition of the single- 
particle overlap Q se \f(t), F(r; t) is replaced by the inverse 
Fourier transform of of the self-intermediate scattering 
function F s (q;t) [H] and Q ra nd is absent. To make con- 
nection with Ref. [5] we chose a = 0.3. 

In Fig. [3]we show the time-dependence of the collective 
and single-particle overlaps for ip = 0.51. For smaller x 
our predictions are qualitatively similar to computer sim- 
ulation results showed in Fig. 2(c) of Ref. [5]. However, 
upon approaching the mode-coupling transition Q(t) pre- 
dicted by the theory shows a classic mode-coupling- 
like two step decay whereas the simulation results ex- 
hibit a continuous increase of both the intermediate-time 
plateau and the long-time plateau. 

In Fig. [4ji we show the length dependence of the so- 
called static overlap, i.e. the nontrivial part of the long 
time plateau of the collective overlap, Qoo — Q ra nd, where 
Qoo = lim^oo Q(t). The length, I = 0.5x _1 / 3 , is the so- 
called confining length [5]. Perhaps fortuitously, for the 
range of lengths corresponding to non-glassy states, the 
values of Qoo — Q ra nd are quite close to those obtained 
from computer simulations. However, according to our 
mode-coupling approach an increase of Qoo — Qrand is 
followed by a discontinuity at the mode-coupling transi- 
tion line whereas the values obtained from simulations 
increase continuously with decreasing I. Importantly, 
growing values of Qoo — Qrand are often associated with 
growing amorphous order and more specifically growing 
static point-to-set correlations. As indicated in Fig. [4^, 
(and also noted in Ref. [5]), for the range of lengths cor- 
responding to non-glassy states, Qoo — Qrand is dominated 
by the small x contribution that originates from expres- 
sion Q . This contribution has a rather trivial origin and 
it should not be associated with any many-body correla- 
tions [3D]. It would be interesting to investigate whether 
there is a similar simple "baseline" contribution for other 
geometries considered in Ref. [5J. Finally, in Fig. [4j) we 
show the length dependence of the relaxation time r de- 
fined through the single particle overlap, Q S cu{t) = e . 
We find that our mode-coupling approach overestimates 
the influence of the pinning on the relaxation time: the 
values of \o{t/t 00 '), where Too is the relaxation time of 
the I = oo (i.e. un-pinned, x = 0) system, are consis- 
tently above those obtained from computer simulations. 
We note that at large I, the length of the relaxation time 
seems to become consistent with the value of the dynamic 
correlation length obtained from the so-called inhomoge- 
neous mode-coupling theory (IMCT) [3"Tl I3"2"] . 
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FIG. 3: Time dependence of the collective overlap Q(t) (solid 
lines) and the single particle overlap Q se u(t) (dashed lines) for 
tp — 0.51 and values of x indicated in the figure. 
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FIG. 4: (a) Dependence of the overlap on the confin- 
ing length I = 0.5i -1 ^ 3 for <p = 0.51. Solid line indicates 
the overlap originating from the small x limit of F(q), Eq. 
Q. (b) Dependence of the single-particle overlap relaxation 
time on the confining length I = 0.5a: -1 / 3 . Solid line shows 
the exponential dependence rn(T/Too) oc exp(— Z/£imct) with 
£imct = 1.66 for ip = 0.51 [32]. Dashed vertical line through 
(a) and (b) indicates the mode-coupling transition. 



In summary, we proposed an alternative mode- 
coupling approach for partially pinned systems. Our ap- 
proach agrees qualitatively with the mean-field analysis 
of the p spin model and with the results obtained from 
the RFOT theory. We showed that the contribution to 
the long-time limit of the collective overlap originating 
from the presence of pinning is dominated by the trivial 
small x limit. This small x contribution is not related to 
any growing point-to-set correlations. 
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